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We investigate the entropy counting for black hole horizons in loop quantum gravity (LQG). 
£N| ■ We argue that the space of 3d closed polyhedra is the classical counterpart of the space of SU(2) 

intertwiners at the quantum level. Then computing the entropy for the boundary horizon amounts 
' to calculating the density of polyhedra or the number of intertwiners at fixed total area. Following 

CN ! the previous work arXiv:1011.5628 [15| we dub these the classical and quantum polymer models for 

isolated horizons in LQG. We provide exact micro-canonical calculations for both models and we 
show that the classical counting of polyhedra accounts for most of the features of the intertwiner 
counting (leading order entropy and log-correction), thus providing us with a simpler model to 
further investigate correlations and dynamics. To illustrate this, we also produce an exact formula 
iy-^ ' for the dimension of the intertwiner space as a density of "almost-closed polyhedra" . 

(N ' 

O" 1 I. INTRODUCTION 

U 

Black holes and their area-entropy law play an essential role in the research of a consistent theory of quantum 
gravity, both as a basic playground to test each model and as a guidance towards understanding the degrees of 
freedom and the structure of the quantum geometry of space-time. Indeed, on the one hand, it is necessary that any 
^ ■ theory of quantum gravity to recover and explain the entropy, radiation and thermodynamics of black holes, while 
on the other hand black hole entropy gives us hints that geometry and gravity could be entirely reformulate in terms 
of quantum information flow. In the context of loop quantum gravity, black hole horizons are described classically as 
isolated horizon [l|. Then general relativity, expressed in terms of triad and connection variables, induces a Chern- 
Simons theory on this boundary horizon. First described as a U(l) Chern-Simons theory in the original work @], the 
analysis was refined and the isolated horizon has been understood to be described by a SU(2) Chern-Simons theory 
[H (see also Q for a review of the formalism and of the relation between the U(l) Chern-Simons theory and the SU(2) 
Chern-Simons theory) . 

At the quantum level, in loop quantum gravity, the state of the geometry of the region outside the black hole 
is given by a spin network state. This state punctures the horizon and induces topological defects on the surface. 
These puncture carry SU(2) representations (spins), which define quanta of area. Then the quantum states of the 
horizon are then given as states of the Chern-Simons theory, or equivalently as the related conformal blocks. At the 
mathematical level, this effectively simplifies to describing the horizon states as SU(2) intertwiner states between the 
SU(2) representations carried by the punctures. Thus, in this context, computing the black hole entropy amounts 
to counting the number of SU(2) intertwiners for fixed total area. In this present work, we will discuss neither the 
conceptual issues of this approach, nor details of the Chern-Simons theory, but we will focus of the entropy counting. 

There has already been substantial work on this topic, calculating the number of intertwiners for fixed values of the 
spins [a. |6| and performing full entropy computations summing over the values of the individual spins for fixed total 
area [3|, l7H10l| using either exact intertwiner counting or canonical methods. Here, we would like to present an exact 
intertwiner counting and focus on providing it with an explicit geometrical interpretation. 

Based on both the work attempting to understand the deeper structure of the intertwiner space [1, [H[ and research 
on defining proper coherent intertwiner states [HI, 0]| , a clear picture of intertwiners as quantized polyhedra appeared 
[il^flil ]. We introduce a classical polyhedra counting and we will focus on comparing it to the true quantum intertwiner 
counting. Our goal is to see how many features of the exact entropy counting can be recovered from a much simpler 
classical counting. In particular, we will be led to an interesting formula of the number of intertwiners for fixed total 
area as an integral over almost-polyhedral configurations. This approach allows to see how one deforms from the 
classical density of states to the number of quantum states. It will appear that the classical and quantum countings 



in 
o 



X 



'Electronic address: |etera. liv ine@cns-lyon.fr 
^Electronic address: danicl.tcrnoO mq.cdu.au 



2 



are very similar and for that most purposes the space of classical polyhcdra can be a perfect simpler substitute for 
the space of intcrtwiners. 

More precisely, we describe the space of polyhcdra with A" faces through the normal vectors to the faces, V\,..,Vn 
in R 3 . The norm of these vectors give the area of the faces while their direction define the planes in which the faces 
lay. These vectors are required to satisfy the closure constraint, V{ = 0. This descri ptio n of classical polyhedra is 
particularly suited to building semi-classical intertwiner states in loop quantum gravity [12Hl4j | . It is well-known that 
such a set of 3- vectors determine a unique (convex) polyhedra, and vice- versa (see 0] for details on the reconstruction 
process). As noticed in (l5j . this can be seen as a polymer chain, made of A" articulated links. This analogy was argued 
to be especially relevant when discussing the energy or angular momentum of the state. Here we will not investigate 
those issues, but we will keep the name and thus refer to counting classical polyhedron as the classical polymer model 
for black holes in loop quantum gravity and to counting intcrtwiners as its quantum counterpart. We will not use 
the model introduced in [l5[ considering classical polyhedra with discrete face area, which can nevertheless be seen 
as half-way between our purely classical and entirely quantum models. 

In the first section, we define precisely our classical polymer model and we describe how to count classical polyhcdra 
for fixed boundary area A and for fixed number of faces N. We perform explicitly the sum over N with a chemical 
potential and we show that the resulting entropy grows at leading order linearly with the area A as expected. We 
discuss the log-correction to the area-entropy law and show that it is directly related to the closure constraint. We 
provide explicit closed formula for both the case with and without closure constraint. We show how to compute the 
fluctuations of the areas of the faces and the correlations between them. We finally show that the number of faces N 
grows in average linearly with the total area A and that this feature is crucial in recovering the area-entropy law. 

The second section reviews and improves the exact intertwiner counting introduced in [8(. In particular, we show 
that we recover the classical polyhcdra counting as leading order at large area and fixed number of faces. Furthermore, 
we obtain the same asymptotic formula for the entropy after summing over A^ up to a shift in the chemical potential. 

The third section investigates the relation between the classical and quantum countings and we provide a formula 
for the dimension of the intertwiner space as counting almost-closed polyhedral configurations, thus reinforcing the 
interpretation of intcrtwiners as quantized polyhedra. 

Finally, the fourth section is a discussion on the role of the number of faces and the chemical potential. Our main 
point is that a determination the chemical potential or equivalently the typical ratio between the number of faces 
and the area will come from an analysis of the coupling of the polyhcdra/intertwiners with the exterior geometry. It 
is the dynamics of our "horizon" with the geometry outside which will tell us that we are not simply looking at an 
arbitrary surface but considering an actual black hole horizon. 



II. THE CLASSICAL POLYMER MODEL 
A. The Micro-Canonical Ensemble and its Phase Space Formulation 

The ensemble that we are interested in is the set of convex polyhedra with N faces and fixed total area A. A 
convenient description for our purpose is through the normal vectors Vi £ M 3 to the faces, where i = 1..N labels the 
face: the vector V, is perpendicular to the plane of the face and its norm gives twice the area \V,\ of the face. The 
total area is then A — i J^. \Vi\. The only constraint that these normal vectors have to satisfy is a closure constraint: 

i 

Such set of vectors determines uniquely a convex polyhedron (embedded in R 3 ). The interested reader will find more 



details in [14J with Minkowski's uniqueness theorem and Lasserre reconstruction algorithm. In the present work, 
we will not describe the reconstruction process of the polyhedron from the normal vector data. We will focus on 
evaluating the number of such convex polyhedra and use it to define the black hole entropy in the context of loop 
quantum gravity. 

We call this the classical polymer model for black holes, following the previous work Indeed, instead of 

interpreting the Vi's as the vectors normals to the polyhedron's faces, one can see them as a linear sequence of vectors 
forming a loop coming back to its starting point (due to the closure constraint) and of given total length (fixed total 
area), thus forming some kind of polymer chain formed by A~ elements. 

Before moving on, we would like to describe the phase space formulation of this ensemble of closed polyhedra with N 
faces. This will become essential when quantizing polyhcdra to SU(2) intcrtwiners and describing the micro-canonical 
ensemble for black holes in loop quantum gravity. 
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Following |16l - [l8| . we introduce ./V spinor variables Zi G C 2 , which are complex vectors living in the fundamental 
representation of SU(2): 



G C 2 



Each such spinor determines a unique 3-vector V(z) G IR 3 by projecting the 2x2 matrix \z)(z\ on the Pauli matrices 



V^{z\a\z), \V\ = {z\z) \z){z\ = \(\V\l + V-a) . (2) 

Conversely, this 3-vector V determines the original spinor z up to a global U(l) phase. Indeed, the spinors e 10 \z) all 
give the same vector, or equivalently said V{z) is U(l)-invariant. 

The next step is to endow the space of spinors with a canonical Poisson bracket [T| 



{zf, zf} = -iSi^, {zf, zf} = {zf, zf} = . (3) 



It is straightforward to check that the components of each of the 3-vectors Vi form (decoupled) su(2) Lie algebra: 



{Vf,VH=26 ij e abc V°. (4) 



Actually the three components V? of the vector Vi are the generators of the SU(2) transformations on the spinor Zj. 
This is the reason why the quantization of this spinorial phase space will be described in terms of recoupling of SU(2) 
representations and that the space of polyhedra will be quantized into the space of SU(2) intertwiners. 
Then the closure constraint J2i V = can be translated in terms of spinors [l3| : 

x;3=o ^n)(^i=^n^)i. (5) 

i i i 

The key remark is that these are first class constraints, generating simultaneous SU(2) transformations on all the 
spinors Zi. 

Finally, the fixed total area condition, \Vi\ = 2A, can also be translated into spinors as ^2Azi\zi) = A, and 
generates global U(l) transformations on all the spinors, i.c simultaneous multiplication by a phase e . This constraint 
clearly computes with the closure constraints. 

At the end of the day, the spinorial phase space (C 2 ) JV //SU(2) x U(l) obtained by symplectic reduction from the 
space of N spinors is exactly the space of (convex) framed polyhedra with N faces and fixed total area up to global 
3d rotations (and global phase transformations). By symplectic reduction, we mean that we impose the constraints 
(closure and fixed total area) and that we quotient by the gauge transformations that they generate. By framed 
polyhedron, we mean that we add the data of a U(l)-phase to each face, which is interpreted as a choice of 2d frame 
on for each face. 

For example, for N = 4, we are looking at the phase space of framed tctrahedra. An arbitrary tetrahedron is 
determined up to 3d rotations by 6 parameters, minus 1 since we are fixing the total area. Then we need to 4 phases, 
one for each phase, minus 1 since we are quoticnting by global phase transformations. In total, this gives 6 — 1+4—1 = 8 
parameters, which reproduces as expected the dimension of the reduced phase space (C 2 ) 4 / /SU(2) x U(l). 

The last main feature of the phase space formalism using spinor variables is the existence of an action of the unitary 
group U(iV) on the space of framed polyhedra with N faces, which leads to a U(iV)-action on the space of intertwiners 
at the quantum level 0, HU, Hi! ■ The action is the straightforward U(A)-action as N x N matrices on the spinors: 

l*> — ► X)0ikl**>- (6) 

k 

This action commutes with the closure constraint, leaves the total area invariant and allows to generate any closed 
polyhedron configuration from a totally squeezed configuration with only two non-trivial faces (with non-vanishing 
area). Let us point out that this action requires the spinor variables Zi and does not exist if considering solely the 
normal 3-vectors Vi- Although the existence of this U(A)-action is very interesting and seems fundamental when 
investigating the structure and symmetries of the space of intertwiners as explained in 0, [Til [T3T ] , we will not need it 
in the present work. 

This concludes this short review of the phase space formulation of the space of convex polyhedra. We will not 
explicitly use the spinorial formalism in the present section, in which we will mainly use the vectors Vi. The spinor 
variables will become relevant when looking at coherent states for quantum polyhedra and the quantum polymer 
model for black holes in section [TVl 
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B. Counting Polyhedra at Fixed Total Area: Entropy and Asymptotics 

Let us now count the number (or more exactly density) of polyhedra with N faces and fixed total area J^i \ Vi I — 2^4: 



pn[A ] = lim / n * (j> - A ^ fe A ■ ( 7 ) 



4ttV5 

where we both impose the closure of the polyhedra and fix the total area. We use the obvious notation Vk = \Vk\- 
The pre-factor 87r is a normalization to insure tat pN=2 = 1, &s we will check later. Let us notice that pn=i = 
by definition. The measure cPV/AirV for each vector is the measure inherited from the analysis of the spinorial 
phase space and SU(2) intertwiners [l6L Il7j|. The parameter e > is a regulator 1 that will be sent to at the end 
of the calculation. It allows to control the behavior of the integrals when the norm of the vectors is sent to oo. In 
principle, this does not matter since the fixed total area condition ^2 k Vk = 2A automatically implies that the norm 
of each vector is smaller than 2A. Nevertheless, we will Fourier transform the constraints and commute the integrals 
to perform the calculation and the regulator will become necessary. As we will see below this regulator plays the 
same role as the +ie shift of the mass in the Feynman propagator in quantum field theory. 
Let us Fourier transform both constraints: 



A' 



-e-Vi 



d 3 Vi f d 3 u . eT , ft f dq 



P 1 J J \ L 4*Vi J (2tt) 3 J 2tt 

= 87r / (0 I e ~ 2iqAl ^ ^ With 7 <* "W e- eV * iqV ^ ■ («) 

The integral I(q, u) converges due to the regulator e > 0. We can compute it exactly as a function of q and the norm 
\u\ by performing an explicit integration over the vector V. We first integrate over the angular part of V and then 
over its norm: 

I(q u) = [ + °°V 2 — e~< v e^ v [ tLj** = [ + °° Vd y e -,v e nV sin ^ = \ (9) 

This is exactly as the Feynman propagator in quantum field theory where u plays the role of the momentum and q 
the role of the mass. Therefore we know how to compute the next integral over u: 

p N [A] = 8n ( ^e~^ A f + °° 4 ™* du 1 



2tt J (2tt) 3 (u 2 - {q + ie) 2 ) N 

8n fdq e - 2iqA 1 ^r^-S) = fdg e _ 2lq (2N 4)! (£Z^, (10 ) 



K 2tt " 2tt 2 4 r(JV) v ; J R 2ir ' (N — l)l(N — 2)! 2 2N ~ A 

The last step is the integration over q £ K, which amounts to the Fourier transform 2 of the inverse monomial q~ n : 

42N-4 

«M°'- 2 ' V-i)Kiv-2)r 

Finally, one can send the regulator e to zero with the pre-factor becoming trivial, e~ 2tA = 1. 



1 By using a little trick on the Fourier transform of the <5-distribution, 2iv&(x) = J dq exp((iq — e)x) for all values of e 6 C, we can avoid 
having to take the limit e — ¥ and have an exact expression for the density of state pjv [A] for all e S R+ : 

where the factor exp2eA will eventually compensate the factor exp(— 2eA) resulting from the integral over the Lagrange parameter q in 
the final expression Hilt . 

2 We use the explicit Fourier transform formula: 



/ 



dq e~ 2i i A _ 2 2n A 2n ~ 1 e~ 2tA 
2^ (e - iq) n ~ 2(2n - 1)! 
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Let us comment on this result. First, one check that the case N — 2 gives us as wanted the normalization 
P2\A\ = 1- Second, the dependence on the total area A is expected. Indeed, starting with the initial definition (0, 
one can do a change of variable V — > V — V/A by normalizing the normal vectors by the total area A. The change 
of variable produces an obvious overall pre-factor A 2N ~ 4 and leaving us with an integral independent from A. This 
area-independent integral leads to the factor (N — 1)\(N — 2)! depending exclusively on the number of faces N. 

Interpreting this polyhedra ensemble as a model for black holes, we have computed the density of polyhedra with 
N faces and fixed total boundary area A. However, (stationary) black holes are supposed to be described classically 
by a single parameter, the horizon area A (or equivalently the mass). Therefore, we need to deal with the number of 
faces. As it was explained in Q, there are two options. One can decide that N is an extra parameter, either a priori 
independent from A — some kind of "quantum hair" as proposed in Q , or fixed a posteriori in terms of the area in 
the case of a physical black hole — as in fl5| as an example. Alternatively, one attempts to sum over the number of 
faces N, in which case it is necessary to postulate the relative statistical weight for N. 

One simple way to reconcile this two points of view is to sum over the number of faces N with a geometric weight 
a N and define the corresponding generating functional: 

P [a,A] ee p 2 [A} + ap 3 [A]+a 2 p 4 [A} + -- = £ a N ~ 2 p N [A} = £ ^f?^ = • ^ 

N>2 N=0 \ + ' V a 

where I\ is a modified Bessel function of the first kind. 

From a mathematical perspective, the parameter a is a weight associated to each face. More precisely, it can be 
interpreted as a rcscaling of the measure d V — > ad V for each normal vector. Since- this measure is actually only 
defined up to a global factor, a simply accounts for this ambiguity. 

From a physical point of view, the factor a" ee e~^ N controls the optimal number of faces and p, can be interpreted 
as the corresponding chemical potential, as suggested in Q . Or we can interpret a as defining a change of area unit 
with A^fa being the physical area corresponding to the algebraic area A. This interpretation is compatible with the 
fact that the factors a N can be produced by re-defining the integration measure as d 3 V — > ad 3 V as mentioned above. 

We define the entropy as usual as the logarithm of the number of states: 

S[a, A] ee log p[a, A]. (13) 

Then one can use the known asymptotic behavior of the Bessel function, or re-derive it by a stationary point approx- 
imation as we will show below in section Til E| to get the asymptotical behavior of the entropy: 

p[a, A] ~ e 2A ^ f - ^(A^)"* +•••), (14) 

S[a,A] ~ 2Ay/^- ^logA+ ... (15) 

On the one hand, we recover a leading order with the entropy proportional to the area, S oc A. As shown in sections 
IA II and IA 21 where we relax the closure constraints, this leading order is simply due to the density of states for N 
random vectors with a fixed sum of their norm (fixed total area). The parameter a = e - ' 1 , interpreted as the chemical 
potential or a choice of area unit, enters directly this leading order. It allows to change the ratio ^ between the 
entropy and area: the entropy depends at leading order on the way we sum over/count the number of punctures. This 
is similar to the mechanism proposed in Q (and earlier in plj ) to fine-tune the ratio | to ^ without changing the 
Immirzi parameter of loop quantum gravity (defining the fundamental area unit with respect to the Planck scale), 
but by playing with the statistical parameters defining the sum over the number of faces. 

On the other hand, we notice that log-correction comes with the expected — | factor, which seems universal. If we 

were considering N decoupled vectors Vi then we would have no log correction. It is the closure constraint J2i ^ = 
that is responsible for the log-correction to the entropy formula. Actually, we have three closure constraints, one for 
each direction of R 3 . Each of them produces a — ^ correction. We illustrate this interpretation in appendix in the 
sections IA II and IA 21 by computing explicitly the entropy for a model without closure constraint and a model with 
closure constraint in a single direction. For instance, if we remove the closure constraint, the log-correction vanishes, 
and if we impose the closure constraint in a single direction (e.g. z) this leads to a — i log-correction. More generally, 
it seems that having n closure constraints leads to a log-correction with factor — 2 . 

On another level, going to the quantum case and counting the dimension of the Hilbert space instead of integrating 
over classical vectors will affect the leading order, but in a manner which will be under control as we will see below 
in the section UTTl Nevertheless the log-correction will still remain the same. 
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C. Spinor Tools: Entropy from Gaussian Integrals 

We have computed above the density of classical polyhedra with N faces for fixed total area defined as an integral 
over N 3-vectors K £ R 3 satisfying the closure constraint J2i Vi = 0. It is interesting to switch back to the spinor 
variables z, <G C 4 used to define the phase space structure of the space of polyhedra. Indeed the integrals above 
over the 3-vectors now become a simple Gaussian integrals over the spinors. This reformulation will be particularly 
interesting when considering the quantum model, where the intertwiner counting will be also expressed as Gaussian 
integrals in the spinor variables over coherent intertwiner states. 

Taking care of the measure when changing variables, we get the following expression for the density of states: 

Pn[A] = Stt Jn^s(^2(z k \z k ) -2A \ 8<® (y2(z k \a\z k )\ 



= 8,e-/ g / J{l d ^ e -^i^-^ . (16) 

We get the same Gaussian integral over each of the spinor variable Zi, which is easy to compute: 

^ e -{z\{e-iq)-iu-3\z) = 1 = 1 ( 17 ) 

7r 2 det2x2 [(e — — iu ■ <r] u 2 — (q + ie) 2 



n 



This allows to reproduce the same expression resulting from the integration over the 3-vectors Vi : 

M1 a 2eA f d q f rf3 " -2«,A ( 1 \ N 

and then perform the integrations over the Lagrange multipliers q and u as before. 

D. Determining the Size of Faces and Probing the Shape of the Polyhedron 

Beyond the mere counting of polyhedra and the resulting entropy for the (classical) polymer model of black hole 
horizons, one can further investigate the finer properties of the distribution over the space of polyhedra. For instance, 
what is the typical size of each face? What is the typical shape of the polyhedron? 

Singling out one face out of the N faces, we can look at the probability distribution for the norm 3 Vi of the normal 
vector. This can be derived from the formula ([5]) for the density of states by computing the integral I(q, u) for the 
(N — 1) other faces while dropping the integral over Vi for our chosen face. Then one can finish the calculation by 
integrating over u and q. However, performing these integrals involve complicated integrals of Bessel functions of the 
second kind. Therefore, we focus on computing the average size of each face and its standard deviation. 

Starting with computing the mean value of the norm Vi of a single normal vector, we insert a factor V in one of 
the integrals I(q,u). This modified integral is easily computed by differentiating I(q,u) with respect to q: 

I(q,u)= [ fL e-^^v^u-v 



AttV u 2 - (q + ie) 2 



T «-o ^ J % - = vX*n ' (18) 

Then we can compute the resulting modified distribution by performing the same integrals as before but with slightly 
different exponents: 

d 3 U dq , iqAt , . ^N-Xf, * „-2eA 2A 2N ~ 3 



p N [A] = 8n ^-^. e -^ A I( q ,u) N -H(q,u) 



(2tt) 3 2tt vy ' ; vy ' ' N!(N~2)\ 



3 Since the measure d?Vi is invariant under 3d rotations, only the norm Vi is relevant while the direction Vi follows the uniform probability 
on the 2-sphere. 
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Then we derive the expected the mean size of a single face: 

<V> = ^} = ^, (19) 
Pn[A\ TV 

which is obvious from the constraint J^. Vi = 2A. Similarly, we extract the variance of the norm V of one normal 
vector using the same technique: 

" / Sr v " ^ w * = = ^TTTW + W^TWf ■ < 20 » 

Replacing I N by I 1 ^^ 1 ! in the calculation for pn[A], we obtain at the end of the day: 



1 ' N{N+l) V[ ' [ > TV V TV + 1 JV»i y/2 [ ' 

which mean that the probability distribution of the size of a face remains fuzzy as the number of faces TV grows. 

We can further look at the correlation (ViVj) between the size of two distinct faces i ^ j. We follow the same steps 
as previously by replacing I N by I N ~ 2 I 2 . This leads to: 

*»• Wi'- ^y+D - (22) 

One can check that, by combining this correlation to the variance (V^ 2 ), we recover the exact relation ^ j{ViVj) = 4A 2 , 
as expected due to the constraint J^i Vi = 2A Furthermore, we can check that the size of the faces become more and 
more decoupled as the number of faces TV grows: 



{VlVj) [V) TV 2 2(TV 2 — 1) TV ^ (V) ^oo U - {Z6) 

By further analyzing the correlations between components of the normal vectors, e.g. (V i a Vj), one would start 
probing the shape of the polyhedron, but this goes beyond the scope of the present work and we postpone it for future 
investigation. 



E. Determining the Number of Faces 

Let us now investigate the probability distribution of the number of faces TV itself. Our density of states is given 
by the series J2 N a N pn[A]. Let us thus analyze the behavior of a N pn[A] in terms of TV and see if it has any clear 
maxima. To this purpose, we use Stirling formula to analyze its behavior at large TV: 

N 2 N 2 ~ ~ 
r A ] JV 2N A 2N-4 AT -2N N r ji Jv 2JV(l-log N)+2N log A a a I / 0/) \ 

PN ^ jv»i 2^ ' PN[A\ ~ ^44 e > A = A^a. (24) 

The exponent ip(N) = 2TV(1 — log TV) + 2TVlogA has a unique fixed point, which is straightforward to determine: 

= d N <p = 2 log A - 2 log TV & N max = A = Ay/a. (25) 

Thus we are in linear regime where the number of punctures TV grows linearly with the total area A. This is consistent 
with the recent proposals [l3j , (l5| and Q , and also with older works on black hole entropy in loop quantum gravity 
(e.g. (H). 

Computing the Hessian at the fixed point d 2 ip\ N = —-«• shows that it is a maxima and that the density a N pjy[A] 

behaves in a first order approximation as a Gaussian around N max = A. Wc check numerically the validity of this 
Gaussian approximation for various values of the parameter. Here figure [T] shows the example of a density peaked 
around N max ~ 20. Using this Gaussian approximation, we can compute the asymptotics of the the sum over the 
number of punctures ^ N a N ~ 2 pn\A\. 

1 / A A 2 p 2 ^ /—= 1 p 2 ^ 

r ,°- [A] . ^/ 2 »_ . ^VTi - w ^ , ,26) 
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FIG. 1: Plots of the true density a N pn[A] for Ay/a = 20 and its Gaussian approximation for N max = Ay / 7x + 1 — 21 (the true 
density is above the Gaussian approximation). The +1 shift comes from the next-to-leading corrections, which can be easily 
computed taking into account the N 2 pre-factor in the exact density. 

which fits perfectly with the asymptotics (fl"4|) for p[a, A] ad S[a, A}. Let us point out that if we had only considered 
the evaluation at the optimal number of punctures PN ma ^[A], we would have retrieved the correct leading order for 
the entropy S[a, A] but not the correct log-correction. 

Finally, we would like to point out the interesting fact that the entropy S[a, A] is simply given by the optimal 
number of punctures N max (up to a constant factor 2) at leading order. Although intriguing at first, this is totally 
expected from the perspective that the entropy is simply counting the number of degrees of freedom, which is directly 
proportional to the number of fundamental blocks thus faces at leading order. 

F. The Linear Regime N — a A 

Instead of summing over the number of faces TV, we could take the point of view that it could be fixed by exterior 
parameters. As we discuss later in section |Vj the number of faces N could be fixed dynamically through the coupling 
of the surface to the external geometry, which would determine if it is a black hole event horizon or an isolated horizon 
or another type of boundary surface. This scenario was already suggested in Q, where the linear regime N = a A 
with a number of faces increasing linearly with the total boundary area was considered. 

Let us start with the density of states p^ [A] and define the entropy sjy [A] = log p^ [A] . A first tentative regime is 
defined at fixed number of faces N. The asymptotics is obvious to write down: 

s N [A] = (2iV-4)logA-log(iV-l)!(iV-2)!. (27) 

We do not even get a leading order scaling proportionally to the total area A. This possibility does not seem to model 
anything close to a black hole or an isolated horizon. 

Following the behavior obtained above for the probability distribution of the number of faces, it seems reasonable 
to consider instead a linear regime where the number of faces N is determined by the total area A and grows linearly 
with it. Let us thus set N = a A with a fixed parameter a. It is fairly easy to extract the asymptotics of the resulting 
entropy using the Stirling formula for the factorials: 

p N= aA[A] ~ 7T-7? (-) 2(TA — > 8*[A] = \ogp N=aA [A] ~ 2<t(1 - log a) A - 2 log A + . . . (28) 

We do recover the wanted area-entropy law with the leading order ratio s/A depending explicitly on the ratio a = N/A 
between the number of faces and the area. A suggestion of Q is to use this freedom in the linear regime to fine-tune 
the ratio s/A to the appropriate physical value by playing with a and without having to fix the Immirzi parameter of 
loop quantum gravity (choice of physical area unit) to a specific value. In this scenario, we hope that the dynamics 
induced by the coupling of our surface with the external geometry will fix the ratio a to some appropriate value. 

A slightly negative point is that the log-correction does not come with the expected — | factor. Nevertheless, this 
issue is easily solved by allowing Gaussian fluctuations of the number of faces N around the optimal value a A. This 
is exactly what we have seen earlier when summing over the number of faces with the "trivial" weight a N defined in 
terms of the chemical potential. We then recover the expected — | log-correction. 
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Finally, if we look at a regime where N grows much faster than A, we lose the linear behavior of the entropy. For 
instance, considering a power law growth N ~ A s with s > 1, the factorials TV! 2 at the denominator of pn [A] will 
dominate over the numerator growing as A 2N and the density of states will decrease as the area grows. 

Thus the only viable regime to recover the area-entropy law is the linear regime where the number of faces N 
grows linearly with the surface area A. This supports considering the sum over N with the simple weight a N , which 
reproduces this exact regime at leading order and indeed predicts a probability distribution for N peaked around an 
optimal value growing linearly with the area A. 



III. THE QUANTUM POLYMER MODEL: COUNTING INTERTWINERS 

A. Computing the Dimension of the Intertwiner Space and Generating Functionals 

At the quantum level, the black hole entropy in Loop Quantum Gravity is given by counting the number of 
SU(2) intertwiners for a fixed total area Q. Considering intertwiners with N legs, we put one spin ji on each leg 
labeled by i = 1..N and we look at the SU(2)-invariant states in the tensor product of the corresponding irreducible 
representations: 

Hf u .. dN = Inv su(2) V^ ® .. ® V™ , (29) 

where V- 7 stands for the (2j + l)-dimensional Hilbert space carrying the irreducible representation of SU(2) of spin j. 

As was shown in [lj, [l6| , and then reviewed in [l4|, G3 > this space of intertwiners comes as the quantization of 
the space of classical polyhedra as defined in section III A[ more precisely it is the space of holomorphic L 2 functions 
on the spinorial phase space. From this point of view, we interpret intertwiners as quantized polyhedra and we see 
the LQG black hole entropy counting of intertwiners as the quantum counterpart of the polymer model defined and 
analyzed in the previous section. For more details on the equivalence between the intertwiner space and quantized 
polyhedra, the interested reader can refer to [1, [U, HH, [I?} and to [3] . 

We would also like to cite [l5j . which proposes an intermediate model between counting classical polyhedra and 
quantum intertwiners. It quantizes the area of the faces while keeping the direction of their normal vector freely 
distributed on the classical 2-sphere. Although this model is not directly relevant to our present work, we give its 
basic definition and a rough analysis in appendix [UJ 

Next, we define the Hilbert space of SU(2) intertwiners for a fixed total area, summing over spin labels ji while 
keeping their sum fixed: 

«J = K,.j»= kv su(2) V*®..®V^. (30) 

Let us point out that we are using here the equidistant area spectrum of Loop Quantum Gravity, where the area 
of a face carrying a spin j is given actually by j in Planck unit (times the Immirzi parameter). This is in contrast 
with the more used spectrum given by the square-root of the SU(2) Casimir operator, \/ j(j + 1) in Planck unit. Let 
us remark that these two spectra only differ by operator-ordering ambiguities (see e.g. j23|). In our context, the 
equidistant spectrum Aj ex j comes naturally from the quantization of the spinorial phase space leading to the space 
of intertwiners 4 . In particular, it allows to realize that the intertwiner space at fixed total area J = ji carries 
an irreducible representation of the unitary group UYA0, which is then interpreted as the group of gauge-invariant 
deformations of N- legged intertwiners at fixed area |8|, LL3| ■ 

The dimension of the space of intertwiners Ti'j was computed in Q actually using the fact that it carries an 
irreducible representation of U(iV) with a two-line Young tableau: 

MJ^^H^-^-f"^- 1 )^^- 2 ), (31) 



4 This equidistant spectrum also turns out to be useful and relevant in other contexts. For instance, in the construction of the EPRL-FK 
spm foam models with Immirzi parameter, one needs to use the equidistant spectrum to solve exactly the diagonal simplicity constraints 
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in terms of binomial coefficients. As also shown in |j§], this dimension can also be computed as the integral over SU(2) 
of the product of suitable characters: 

N 

d N [J] = dg £ ~[xj t (9), (32) 

SU(2) Y.tii=J 1 

where dg is the normalized Haar measure and Xj the SU(2)-character of spin j (i.e the trace of the matrix representing 
the group element g in the representation of spin j). 

For a fixed number of faces N and a large area J — > oo we get 

j2N-4 jyj2N-5 

dN[J] J^oo (JV-l)!(JV-2)! + (N-l)\(N-3)\ + ' ' " (33) 

which in the leading order is the same as the classical density of states pn[J] derived in (fTTj) . Thus the number of 
states djy[J] and the density of states pjvLA] at fixed number of faces N coincide for large areas is consistent with 
seeing the intertwiner counting as the quantized version of the classical polymer black hole model counting classical 
polyhedra. This is a strong result showing the consistency of the approach of counting classical polyhedra for modeling 
the isolated horizon state counting in loop quantum gravity at least at leading order. As a consequence, we expect 
very similar results for the generating functional and the resulting entropy. 

The next step is to sum over the number of legs N (or equivalently the number). At fixed total area J, we sum 
over N with a geometric weight a and obtain a hypcrgcomctric function: 

d[a,J]=^^d N [J] = E^ ^'^fJ+lV = + + *,*,*)■ (34) 

N>2 N=0 ' V 

As shown in Q, or using the known asymptotics for the hypergeometric functions, or working out a saddle point 
approximation for the above series, one can extract the large area behavior of this generating functional: 

S[a, J] = log d[a, J] j ~ i -2Jlog(l - yfa) - ^ log J. (35) 

Thus we find that classical polyhedra and quantum intertwiners lead to very similar entropies, which differ only in a 
renormalization of the weight a associated to each puncture (or equivalently the chemical potential): 

S[a,A] ~ lAyfa- ^\ogA — > S[a, J] ~ — 2 Jlog(l — y/a) — ^ log J. (36) 

Identifying the total area A = J with the sum of the spins, we see that the ratio S/A differs between the classical 
model and the exact quantum computation. This change of leading order between the classical polymer and quantum 
polymer is due to the quantization and subsequent discretization of the area spectrum. Both entropies have the same 
leading order when the weight associated to each puncture is very small, a — > 0. Moreover, simply renormalizing the 
chemical potential by changing \fa into — log(l — \/a), allows to get the quantum entropy from the classical entropy 
formula. We will this classical <-» quantum transition later in Sec. IIV1 Let us also point out that the quantum entropy 
S[a, J] is always smaller than its classical counterpart S[a, J], 

Moreover an interesting point is that the log-correction is exactly the same in both cases. This robustness of the log- 
correction can be interpreted as coming from the number of constraints that we are imposing: the closure constraints 
in R 3 in the classical case corresponds to the requirement of SU(2)-invariance in the quantum theory. Since SU(2) 
is also 3-dimensional, the SU(2)-invariance translates into 3-dimensional constraints, which lead to the factor — |. 
Wc illustrate this by working out explicitly in appendix [B] the entropy counting for models without SU(2)-invariance 
(where there will be no log-correction) and with a simpler U(l)-invariancc (where we recover a — ^ log-correction). 

Finally, using the Stirling approximation for the factorials for large J and N, we can compute the optimal number of 
faces N max given fixed total area J and chemical potential a. Identifying the maximum of the probability distribution 
a N dN[J] easily gives at leading order: 

N „ J ^ ( 37 ) 

1 - Ja 



which is almost the same as the classical maximum up to the factor 1 — y/a in the denominator (which nevertheless 
becomes trivial if a is sent to 0). 
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B. Computing the Generating Functional 

Let us look into defining the generating functionals for the dimensions d^[J]. This will be particularly useful when 
comparing the classical and quantum entropies using coherent states in the next section. Following @ , we can define 
a generating functional /jv[i] = X^jgN t 2J d-N[J] by summing over the total area J. We can express it as an integral 
over SU(2) using the formula (j3"2")l for dp/ [J]: 



IN 



[J}= [ dg J2 UxiM fN[t] = E^W] = / dg I £ t^(g) 

su(2) £f jw * JeN (2) V eN / 2 



AT 



(38) 



Taking into account that the characters are class functions (invariant under conjugation) and parameterizing the 
group clement g in terms of its class angle (f>, we write the generating functionals /jv[i] as integrals over the single 
angle (p: 

t m 2 r • 2 ... ^ ,2, s in (2j + 1)<A 2 /■"• sin 2 

^--y^'w^^-s^j = -| # ^_^- TF . (39) 

Then we can introduce an overall generating functional by also summing over the number of faces N: 

F[at,t] = a N f N [t]=J2 t2J £ =^( 2J (l + a 2 4«> J]), (40) 

W6N ./ NeN J 

which is thus also a generating functional for the dimensions d[a, J]. Explicitly computing the sum over N allows to 
write F[a,t] as a trigonometric integral: 

r n 2 r 2 , , , 2 - 2iCOS0+ 1) , , 

F[a,t} = - sm 2 4>d4>-i — , Zs ■ 41 ) 

Trio (t 2 ~2tcoscj}+l)-a v ' 

This integral can be computed explicitly from a,t ~ Q: 

F(a,t) = 1 + | + ^(l-a-V(«-(l+*) 2 )(a-(l-t) 2 )) (42) 

For < a < 1, this generating functional is well-defined for t around 0. Then the lowest (non- vanishing) pole is at 
t c = 1 — s/a, or equivalcntly at a = (1 — t c ) 2 . This gives us the asymptotic behavior for the coefficients of the power 
series in t 2J 



log a N d N [J] - -2Jlogt G = -2Jlog(l - Va), (43) 
which reproduces exactly the leading order of the asymptotics (|35p of the entropy. 



C. A Remark on the "No-Trivial Puncture" Model and Entropy 

An issue with the intertwiner counting presented above in section A[) is that it includes vanishing spins on the 
intertwiner legs ji = 0]. This corresponds to faces with vanishing area. Therefore, one would like to ideally remove 
them and perform a "no trivial-spin" counting. As was shown in Q and will be explained below, the "no trivial- 
spin" entropy counting can be re-absorbed in a re-definition of the weight (chemical potential) a. Indeed removing 
trivial spins amounts simply to a specific summing over the number of faces N. Thus all one needs is the generating 
functional d[a, J] summing over the number of punctures TV and this density of states is enough to take into account 
removing 0-spins or even more general restrictions one would like to impose on the intertwiner counting. In particular, 
the "no trivial-spin" entropy Sz[a, J] will have the same asymptotical behavior as S [a, J], with a renormalized factor 
for the leading order and the same log-correction. 

Before showing the details of the computation of this "no trivial-spin" entropy S0 [a, J], we would like to discuss 
the legitimacy of removing the contributions of 0-spins. From a classical and geometrical point of view, it is natural 
to remove faces with vanishing area, since these arc simply not actual faces. From a quantum point of view, removing 
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O-spins amounts to removing the ground state of the Hilbcrt space and the legitimacy of doing so can be questioned. 
First, considered solely the horizon, there is no principle or reason to exclude quantum superposition of spins including 
0-spin on a face/leg of the intertwiner. Second, considering the coupling of the black hole horizon with the exterior 
space-time geometry, the number of faces N is a property of the quantum state of geometry of the external region 
and should come from the edges of this spin network state puncturing the black hole horizon. Nothing forbids these 
edges to carry superpositions of various spins including 0-spin. From this point of view, the number of faces N should 
be dynamically derived from the coupling of the black hole with the external geometry. However, to go further in this 
direction would require a kinematical and dynamical model of the near-horizon geometry in loop quantum geometry 
and it is not enough to have a kinematical picture of the horizon decoupled from the rest of the space-time geometry. 

The number of intertwiner Dn[J] for fixed total area J with a number of faces N constrained to carry non- vanishing 
spin is easily related to the number of intertwincrs djv[</] with possible O-spins through the binomial transform, as 
used in Q: 

dw[J\ = J2 (k) D kW> D K [J]=Y,i-l) K - N (n) MA (44) 

K=0 ^ ' JV=0 ^ ' 

in terms of binomial coefficients. Then using the series expansion 5 of (1 — x)~ N , one can easily relate the two generating 
functionals: 

u N d N [J] A aN ° N t J ] 

JV JV 



oo oo oo , s 1 , s JV 

^a«D K[ J] = £ (-!)*""«* (n)= T^E(t^) MJ]. (45) 

K=0 JV=0 K=N V 7 JV V 7 

Therefore, summing over faces carrying non-trivial spins is fully equivalent to summing over faces with arbitrary spins, 
up to a controlled renormalization of the weight a of individual faces. 

Thus, introducing the non-trivial spin entropy S [a, J] = log J2n>2 aN ~ 2 ^N[J], where we removed the nonsensical 
case N = and the vanishing term D\ [J] = 0, its asymptotics at large area J ^ 1 is automatically derived from the 
known behavior of the entropy S [a, J] given by ([351 : 

5 [a ) J] j ~ i -2Jlog(l-y^)-flo g J + .... (46) 

To conclude with the issue of removing trivial faces or not, we do not need to consider explicitly the case of removing 
the contribution of O-spins since this is a possibility already taken in account by considering the generic generating 
functional d[a, J] summing over the number of faces N. 



IV. BRIDGING BETWEEN THE CLASSICAL AND QUANTUM MODELS 

The Hilbert space of SU(2) intertwiners with N legs is understood to be the quantization of the space of closed 
polyhcdra. It would thus be enlightening to understand how the intertwiner counting can be actually interpreted as 
computing a number of classical polyhedra or "almost-classical" polyhedral configurations. To this purpose, we will 
use coherent intertwiner states to write a semi-classical formula for the dimensions of the intertwiner spaces as an 
integral over fuzzy polyhedra. This would allow to bridge between the entropy formula for the classical and quantum 
polymer models described previously. Before getting into the details of coherent states, let us start with a simple 
remark on the entropy counting. Indeed, as we already pointed out earlier, it appears that we can obtain the classical 
entropy formula as a double scaling limit of the quantum entropy formula sending both the area to infinity and the 



5 We use the following scries expansion around x = 0: 



x _ k 

Cl _ x )N+i ~ x 

{1 X > K>N 
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chemical potential to 0. For fixed parameters A and a, we set J = kA and a = k 2 a in the intertwincr counting and 
send the scaling parameter k to oo: 

S[a,J] ~ -2Jlog(l-Va) — >S[k~ 2 &,kA] ~ -2«Alog(l - n^V^) ~ 2AV& ~ S[a, A] (47) 

Let us point out this double scaling limit is the one that keeps fixed the most probable number of faces N max of the 
intertwiner/polyhedron. Thus in a sense, we rescale the chemical potential at the same time as the area in order to 
keep the number of faces finite. 

This suggests that we should be able to recover the classical polymer model in a certain large area regime of the 
quantum model, or vice-versa we should be able to interpret the quantum polymer model as a suitable deformation 
of the classical model. 



A. Semi-Classical Formula for Intertwiners and Quantum Polyhedra 

Here we would like to propose an alternative expression for the number of intertwiners as an integral over classical 
vectors representing fuzzy or actual polyhedra. Following the approach of [l2l . ll4T | to define coherent intertwiners 
through group averaging, we can introduce the SU(2) coherent states \j, z) labeled by spinors z € C 2 living in the 
fundamental representation of SU(2). Those \j, z) states form a over-complete basis of the irreducible representation 
V 3 of spin j and transform covariantly under SU(2) (see e.g. [HI, [l|| fl7l . [lj| for details on their definition and 
properties). Here, we will use the decomposition on the identity on the representation V 3 of spin j: 



e~^d A z 



\3,z){j,z\, (48) 



where (z\z) stands for the norm of the spinor z in C 2 . This allows to express the characters Xjis) as Gaussian integrals 
over C 2 : 

Xj{9) = m J M9\3,z) = j^y j (z\g\z) 2j , (49) 

where g\z) is the natural action of the group element g on the complex spinor z in the fundamental representation of 
SU(2). We can insert this expression in the definition of the dimension cIn[J] as an integral over SU(2) of a product 
of characters: 



N „ N 



d N [J]= [ d g ]T f[xi<(g) = [f[ e Zi l dZi E rnW I dg{zMzi) 2j >- (50) 

Jsu(2) j2 N ji=J * ' J2 N ji=J 

We can write this integral in term of the 3- vectors Vi = (zi\a\zi) (since the expression is invariant under multiplication 
by a phase of each spinor Zj). Parameterizing explicitly the group element as g = cos 9 + i sin 9u ■ a in term of the 
class angle £ [0,7r] and the rotation axis u G <S 2 , we have: 

d N [J] = j E nS yy- J -n 2 OdO — ( Vi cos6 + i S in9u ■ . (51) 

Besides the sum over the spins ji, this is an entirely classical expression for the dimension of the intertwincr space. 
Following the original approach of [l2j], one can do a saddle point analysis of the integral over SU(2) for large spins 
ji, in which case one finds that this integral is exponentially suppressed except if the closure constraint 'Yl li jiVi = 
is satisfied. This is not the closure constraint of the classical model. We can nevertheless restrict the integration 
of over the coherent states to unit spinors, (zi\zi) = 1, and thus to unit vectors Vi G S2 which we write with a 
small v to distinguish them from the normal vectors Vi with arbitrary norm. These still provide a decomposition of 
the identity on the SU(2) representation spaces V 3 . This are actually the procedure used in [H, [13] when defining 
coherent intertwincr states. Then the closure constraint determining the saddle point at large spins is ^^jiVi = 0, 

which can be written simply as Vi = after defining the normal vectors as V, = jiii with Vi = ji and thus 
Si = 12iJi = J ■ The problem with this formula is two-fold. First, the norm of the normal vectors are quantized 
and discrete and we do not have an integral over classical continuous vectors in R 3 . This corresponds more to the 
polymer model by Bianchi a la Bohr-Sommerfeld as presented above in section [C] Second, the integrand does not 
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admit a simple natuural expression in terms of these normal vectors Vi with discrete norm equal to j, . Instead, we 
would like to propose here another semi-classical formula for the intcrtwiner counting based on improved coherent 
intertwiner states, as defined in [H, H3: which are constructed as U(7V) coherent states. 

Indeed, it was proven in Q that the Hilbert space of intertwiners with N legs and fixed total spin J = ^ j, is an 
irreducible representation of the unitary group U(iV): 

U N i = Inv su(2) V™ ® ..®V> N . 

It was then proposed in [l3[ to define coherent intertwiner states in T-L jj consistent with the U(iV) action. Here we do 
not need the V(N) action, we will thus refer the interested reader to [lj, [I?} for details on the definition and analysis 
of those U(iV) coherent states. We will simply use the formula for the decomposition of the identity on the Hilbert 
space % J N : 

4 = jy^jVl)! /II^? 6 " 1 ^ UM})(JM}\, (52) 

where the coherent intertwiners are labeled by the total spin J and N spinor variables £ C 2 . Those states have a 
non-trivial norm: 

This provides a beautifully simple formula for the dimension of the intertwiner space: 



2 



(54) 



where we switch the integral over the spinors to an integration over the normal 3-vectors. 

It is interesting to check this formula explicitly. We use similar techniques as when computing the classical density 
of states pjv[^4]- We introduce the Fourier transform parameter u and q and write: 



d M \J] - - fjJ e V,d3Vl (B d -d 2 ) J e^^^^qT.* 



u—q—Q 



where I(q,u) is the same integral as computed in section Hi Bl 



(55) 



J M = / rl e-^e*^ = ^ , (56) 

J 47TK M Z — (<7 + ip) Z 

where we will take j3 = 1. It is easy to compute the action of the Laplacian on the integral kernel: 

(d Ui 8 Ui - flg) 1 w = ~ mN - 1] N+1 . (57) 

By a straightforward recursion, we get: 

f BB j?)J 1 _ lt2J (N + J-l)\(N + J-2)\ 1 

1 " < U< " j (iP_ {q + ij3) 2f - ^ (N-l)\(N-2)\ {U 2_ {q + m 2f+J- W 

Putting all the pieces together and setting back j3 = 1, we get: 

, r T] (N + j-iy.(N + j-2y. 

dN[Jl ~ J\(J+1)\(N-1)\(N-2)V (59) 
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which is exactly the expected result as given in ean. (|31l) . We could have also compute this integral by directly 
performing the Gaussian integral over the spinor variables, in which case it would have turned into a combinatorial 
problem of counting the number of pairings in the polynomial (J, {zi}\ J, {zi}} using the Wick theorem. 

We therefore have a exact formula for the intertwiner counting as an integral over classical 3- vectors Vi. But we 
do have neither closure constraints nor a total area constraint. The natural question is whether this integral can be 
interpreted as an integral over semi-classical configurations or fuzzy polyhedra. 

Looking more closely at the factor ((X)i^) 2 — IX)i^| 2 ) J ; we can factor out the total area J^V* an d we get a 
positive expression smaller or equal to 1 to the power J: 

f, \T,iVi\ 2 Y 



Since the normal of the "closure vector" ^ . Vi is always smaller or equal to the total area J2 i Vi , this expression goes 
to as J grows except if and only if the closure constraint is satisfied Vi = 0. Thus our integral is clearly peaked 
on closed configuration and more so as the total spin J grows. And we do have an integral over some kind of fuzzy 
polyhedra. We can be slightly more precise and write: 

i i i i 

where it becomes clear that we have approximatively at large area J a Gaussian peaked on the closure constraint for 
classical polyhedra. 

On the other hand, the total area Vi does not seem to be fixed. In order to find out, let us compute the average 
of Vi and its fluctuations. We can compute directly the average of e 17 ^ Vi and this will contain the average and 
all the higher moments. This does not require any further calculation since it amounts to taking /? = (1 — 77) and 
plugging it in the expressions above. This leads to: 

(pv^v,) = 1 /"tT— P -0--v)T.iVi 



/n^- (i -^ (E*) a -iE3i 



d N [J]2 2J J\(J + 1)! 

1 _ v? ?" (n + N + J-l)\ 

(1-V) N+J ^n! (7V+J-1)! • [ > 



This perturbativc expansion allows to read the average and the fluctuation: 



(2> =(,+»), v/gggpgsgg 



(62) 



Thus the total area Vi is peaked as the total spin J grows, but its average value has an intriguing +N shift. 

At the end of the day, we do have a semi-classical formula for the intertwiner counting as an integral over the normal 
3-vectors representing some fuzzy polyhedra almost satisfying the closure constraint and almost of fixed total area. 
This provides a simple geometric interpretation to the counting of intcrtwincrs. Using this semi-classical formula, one 
could now go further and compute all the averages and correlations for the normal 3-vectors in the quantum polymer 
model. 



B. Gauge-Fixing to Closed Polyhedra 

Up to now, we have written the dimension of the intertwiner space as an integral over almost-closed configurations 
representing some kind of fuzzy polyhedra. It is actually possible to gauge-fix this integral to exactly closed polyhedral 
configurations. Indeed following the previous work on coherent intertwiners [l3|, 53 j it turns out that the coherent 
intertwiner states are invariant under global SL(2, C) transformations on the spinor labels Zi, while this affects the 
closure constraint: 



\ Zi ) eC 2 ^A|z,), \J,{zi}) = \J,{A Zi }}, \zi)(zi\ ->■ A|z i )(z,|A t , 



(63) 
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where A G SL(2, C) acts as a 2x2 matrix. The key point is that, for an arbitrary set 6 of spinors Zi, it is always possible 
to find a transformation A such that the transformed spinors satisfy the closure constraint, A J^i oc I. 

Indeed, as noticed in fl3l ]. the 2x2 matrix X = ^2n\zi){zi\ is Hermitian and strictly positive and can always be 
written as X = AA _1 A -1 ^ with A > and A G SL(2,C). Reversely, starting with arbitrary closed configuration, 
we can obtain any random (non-closed) configuration by SL(2, C) transformations. This way, we can gauge-fix the 
integral over random configurations to an integral over exactly closed polyhedral configurations (up to a measure zero 
set): 

MJ] = JHj\ IV /n^T* fe !*><*! - ^Efol*) 1 ) (J,{*}\M«}) I e-Ei<"l AtA l"> 

■\ "r )■ J i \ j i / JSL(2,C) 

n£l^(^})^E^(E^, (64) 



2 2J J!(J+1)! J 114^ 
where the measure ^({V;}) reflects the volume of the gauge orbits under the SL(2, C) action 7 : 

**({?<}) = / e-iE^'iVAtA = 

-/SL(2,C) 

where we have assumed that the and Vi satisfy the closure constraint. This measure factor is a function of only 
the total area Vi, we will not compute it but we refer the interested reader to similar calculations in [IH, [2ol [26j. 

Our purpose here is merely to point out the existence of such a formula giving the dimension of the intertwiner space 
as an integral over closed classical polyhedra but with a non-trivial measure. However, the integral over almost-closed 
configuration derived above is simply more convenient to handle for practical purposes. 



C. The Generating Functional as Integral over Fuzzy Polyhedra 

Interesting, we can also give an interesting geometrical interpretation for the generating functional fx [t] defined in 
as the sum over the total spin J. Starting with the expression (|50[) for the dimension of the intertwiner space 
in terms of SU(2) coherent states, the sum over J allows to sum over all the individual spins ji and finally obtain a 
simpler formula: 

f N [t) = ^dAJ] = [ dg [fl^e-We^'M (65) 

Switching the integration from the spinors to the 3- vectors V, and parameterizing the group element g = cos 9 + 
isinOua as before, we get: 

f N [t] = [ TT e ~*ff* f-sin 2 ^/ d 2 ue cose ^e isine ^^. (66) 
J ax 4 n Vi Jo tt J S2 

Written as such, the integration over SU(2) truly resembles imposing the closure constraint J2i Vi> U P to the fact that 
the Lagrange multiplier p = sin Oil has a norm bounded by 1. We can thus interpret this formula as an integration 
over "almost-closed polyhedra" . 



6 The only exception is the degenerate case when all the spinors are proportional to each other. In this case, the 2x2 matrix X = \zi){zi \ 
has a vanishing determinant and there does not exist a SL(2, C) matrix mapping this degenerate set of spinors on a closed configuration. 
However, this set is of measure zero and does not play any important role in our integrals. 

7 Writing explicitly the SL(2, <C)-matrix A in terms of its four matrix elements a, b, c, d we can express the measure as a quasi-Gaussian 
integral in terms of A = i Vi: 

f,({V t }) = J d 2 ad 2 bd 2 cd 2 d8^(ad-bc-l)e~ A ^ 2 +W 2 +^ 2 + ^ = f^f d 2 ad 2 bd 2 cd 2 de^ T -^ ad ' hc ' 1 '> e 'A(\a^ + \^+\^ + \^) _ 

where e is an arbitrary positive regulator (the integral does not depend on its value) and the integral over the a, b, c, d gives simply the 
inverse square-root of the determinant of the Hessian. Then it still remains to compute the integral over T. 



17 



If we were to perform explicitly the integration over SU(2), we would actually get the expression derived using the 
U(iV) coherent states: 



where I\ is a modified Bessel function of the first kind. However this expression does not make clear that the 
integration is peaked on closed polyhedral configurations. 



V. ON THE ROLE OF THE NUMBER OF FACES 



We have discussed a classical polymer model for the black hole horizon, modeling it as a classical polyhedron with 
N faces. Quantizing this into a quantum polymer model leads to the description of quantum isolated horizon whose 
states are iV-valent intertwiners in loop quantum gravity. We have computed the number of states for fixed area, 
giving a density of states /0./v[A] in the classical model and the number of intertwiners cIn[J] where A and J are the 
total boundary area in the corresponding model. Then we would take the logarithm of this number of states and get 
the horizon entropy. An obvious issue comparing the standard classical description of the black hole entropy is that 
we now have two numbers defining the properties of the horizon, the total area A (or J) and the number of faces (or 
legs of the intertwiner) N. We thus need to comment on the role of N. 

In the recent work (9( , it was proposed that N be considered as some "quantum hair" , needed to define the properties 
of the black hole horizon at the quantum level and without any classical counterpart. Then, as was shown in [8|, if we 
want the entropy log/3jv[-<4] or log d^[J] to scale at leading order linearly with the total area, we need the number of 
faces N to grow linearly with the total boundary A. In this case, the entropy/area ratio depends on the N/A ratio. 
One can thus fine-tune this ratio in order to ensure the correct entropy/area ratio. Going further it was proposed in 
3 that the ratio N/A defines various regimes of the theory and should be selected by the actual quantum gravity 
dynamics or by the physical requirement of the surface to be a true black hole horizon. For instance, a black hole 
horizon could have a certain N/A ratio while another type of surface has a different value of this ratio due the different 
quantum gravity dynamics in the neighborhood of the surface. 

In such a scenario, the dynamics of the near-horizon geometry plays an essential role. Indeed the number of faces 
N is actually in loop quantum gravity a property of the spin network states describing the exterior geometry and 
puncturing the horizon surface. It is the near-horizon geometry state and its dynamics that determines that the 
considered surface is actually a isolated horizon. The study of the near-horizon geometry plays actually an essential 
role in other approaches to the black hole entropy, such as the conformal field theory calculations (e.g. [27j ) or the 
more recent proposal [28[. From this perspective, the near- horizon geometry state will determine what are N and 
A. It is not unreasonable to expect that the specific shape and choice of the near-horizon wave-function leads to a 
relation between N and A. Then in order to discuss the value of iV we would need to enlarge our polymer model for 
the horizon to take into account its coupling to the exterior geometry. 

Another possibility is to consider that we have already taken into account and integrated over the dynamics of the 
exterior geometry and that it leads to a certain weight for the creation and annihilation of a face/puncture on the 
horizon. For this point of view, it is thus natural to sum over the number of faces TV together a certain statistical 
weight a N reflecting the probability of the horizon of having TV faces. Such a weight implicitly assumes that the 
creation/annihilation of a single puncture is independent from the other faces. This might have to be revised once we 
have a decent model of the dynamics of the near-horizon geometry. Putting this aside, we are led to an entropy as 
we have defined in the present work S [a, A] = log^jy a N pn[A]. 

The parameter a is the chemical potential and controls the probability distribution of the number of faces. We 
can actually look at this probability distribution in our polymer models and we easily find that for large areas it has 
approximatively a Gaussian shape peaked on a optimal value N max which depends on the chemical potential a and 
grows linearly with the horizon area A. This means that we actually recover effectively the same picture as in the 
previous scenario with the number of faces fluctuating and being determined dynamically to a most probable value 
growing with the area. Then the ratio N max / A depends on the chemical potential a, which reflects the strength of 
the coupling between the horizon and the outside geometry. 

Finally, in order to go further, we will need to go beyond the present kinematical description and look at the 
dynamics of the horizon and its coupling to the outside geometry state. Such a proposal was recently put forward in 
|28| where a quantum Hamiltonian for a quantum Rindlcr horizon in (covariant) loop quantum gravity was defined. 
This should not only lead to a better understanding of the role of the number of faces/punctures but will most 
importantly necessarily lead to a description of the thermal radiation and black hole evaporation process. 
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Conclusions 



Following the identification of black hole boundary states as SU(2) intertwiners and the the description of those 
intertwiners as quantized polyhedra, we introduced a classical polymer model, which defines the black hole entropy as 
a counting of classical polyhedra with fixed total area, and we have explained how its quantization leads to a quantum 
polymer model reproducing the usual LQG intcrtwiner entropy counting. This shows that most of the features of the 
intertwincr counting arc already present in a classical computation of the density of polyhedra at given area. 

First, summing over the number of faces of the polyhedra, or equivalently the number of punctures on the horizon 
or the number of legs of the intertwiners, leads to a leading order entropy scaling with the total boundary area. 
The entropy-area ratio depends on the value of the chemical potential, which also controls the optimal number of 
faces in terms of the boundary area. We argue that the precise value of the chemical potential should be determined 
dynamically through the coupling of the polyhcdra/intertwiner with the exterior geometry. 

Second, we recover the expected — | log-correction. We show that it comes directly from the closure constraint 
imposed on a set of N vectors (in K 3 ) for them to properly form a geometrical (convex) polyhedra (as normal vectors 
to the faces). This constraint is the classical counterpart of the requirement of SU(2) invariance at the quantum level, 
which defines intcrtwiner states. By a full and explicit analysis, we remove and relax the closure constraint and show 
that it affects the log-correction (but obviously not the leading order). We do the same at the quantum level and 
show how removing or relaxing the requirement of SU(2) invariance changes the entropy at next-to-leading order. 

Finally, wc provide a semi-classical formula for the intertwiner counting as a density of almost-closed polyhedral 
configuration, which shows how to interpolate from the classical model to the fully quantum model. 

As next steps, we see two directions. At the kinematical level, we could dig further than the mere entropy which 
counts the total number of states, but we should compute the correlations between the faces (or equivalently the 
punctures on the horizon or the legs of the intertwiners) . This should allow us to probe the shape and structure of 
the horizon, thus probably understand how the angular momentum and higher multipolar moments (see (29| for the 
definition of multipole moments of isolated horizon) can be taken into account as suggested in [l5| . But this could also 
give us some hints to the evaporation process, which seems to be related to the entanglement between the punctures 
on the horizon as suggested in [f|. Then we should go dynamical, model and study the near-horizon geometry and its 
coupling to the boundary horizon, in order to get a clear picture of the thermodynamical role of the entropy which we 
computed and the resulting Hawking radiation. This would truly validate, or invalidate, the physical interpretation 
of the mathematical polymer model for black hole horizons which we have discussed here. 
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Appendix A: Weakening the Closure Constraints for the Classical Polymer Model 



1. Ensemble Without the Closure Constraint 



Let us now investigate the effect of relaxing the closure constraint. Let us thus fix the total area and introduce the 
following density of states: 



with the regulator e as before. We use the same technique as before. We take the Fourier transform of the ^-distribution 
and swap the integrals: 



P%[A} = J ^.e-^ A l\q) N with J°(g) = j 



?S ■ (A2) 

The integral kernel 1° is simpler to compute than in the case with the closure constraints: 

r+oo -, i 

I°(q) = / VdVeT&eW = -—-^ = ? ^ ■ (A3) 
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Wc can then perform the integral over q and get: 

<A4) 

The factor e~ 2eA becomes trivial when the regulator is sent to e4 + . The area dependence A 2N_1 is easily derived 
by doing the change of variable V — s- V/A in the definition (|A1|) of the density of state p°. 

As in the case of closed polyhedra, we define the generating functional p°[a, A] by summing over the number of 
punctures N with an arbitrary geometric weight a": 



DC 



(r, A \2N+1 

p°[a,A] e, J2 <x N Pn[A] = £ ^ N+1 \^ = v^smh(2AV5). (A5) 

N>1 ^=0 ^ '' 

The entropy is defined as the logarithm of the density of states and it is straightforward to extract its asymptotical 
behavior at large area: 

S°[a, A] = log p°[a, A] = log sinh(2A v / a) ~ 2Ay/^ - log 2 . (A6) 

A first remark is that we obviously check that the entropy without the closure constraint is larger than the entropy 
computed earlier with the closure constraints, S°[a, A] > S[a,A]. Then the point which we wish to illustrate is that 
we do not have any log-correction term and that the asymptotical formula goes directly to the constant correction 
term. This is coherent with the interpretation that the log-correction is specifically due to the closure constraints. 



2. Weakening the Closure Constraint: z-Closure Model 

Let us now introduce and analyze a model where we impose the closure constraint along a single axis, J2i = 0- 
This can be considered as the classical counterpart of the "old" calculation for black hole entropy calculation in Loop 
Quantum Gravity that is based on counting U(l) intertwiners 0. 

We define the following density of states at fixed total area A and satisfying the z-closure constraint: 

&w - 1 ft 5 (e % - 2 ^) * (e yt) ( A7 ) 



Taking the Fourier transform of the constraints it becomes 



PM = J ^ ge"** 1 J 1 (?,«)" with l\q,u) ee J ^e-W^" • (A8) 

This new integral I (q, u) is actually equal to the integral I(q,u) with the full closure constraints: 

I'M = fwe-^-^. 2 1 (A9) 
Jo u — {q + it) 

The difference between the expressions (|A8[) and ([5} is that we are now integrating over a single variable ael instead 
of a 3- vector ael 3 . This gives: 



^ - L d i^" A C d i 



i 



(u 2 - (q + ie) 2 ) N 



dq & -2 iqA (e-iq) 1 -™ ( 2N -2 



2tt 2 2JV -! \ N-l 

A2N-2 

= V~*vhw < A10 » 

As before, the factor e~ 2At vanishes as the regulator e is sent to 0. The power of the area, A 2N ~ 2 , can be again 
simply derived by performing the change of variable Vi — > VI /A in the initial definition of the density pjyL'l]. The 
difference with the case with the full closure constraints, besides the slight shift in the power of A, is the change of 
the (N — 2)!-factor with a (N — l)!-factor . This little difference will be responsible for a change in the log-correction. 
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Indeed, we now define the generating functional by summing the number of faces ./V with the geometric weight a . 
p 1 [a,A}=pl[A}+apl[A]+a 2 pl[A} + ... = £ a^pM = - £ a N w = -/ (2^), (All) 

N>1 " JV=0 

where Io(x) is the modified Bessel function of zeroth order. We take the logarithm of this density of states to define 
the entropy. It is then straightforward to derive its asymptotic expansion: 

S 1 [a, A] = log p 1 [a, A] ~ 2 Ay/a - - log A - - log a - log4 - - log -it + . . . (A12) 

We find the same leading order as with and without the full closure constraints. The main difference is with the 
next-to-leading order, where the log-correction now comes with a — ^ factor. 



Appendix B: Weakening the Closure Constraints at the Quantum Level 



1. Without the Closure Constraint 



Let us start with the quantum polymer model and remove the closure constraints as in the classical case (section 
IA lj) and see what happens. Now we do not require the invariance under SU(2) and our Hilbert space for N faces and 
fixed total area J is: 

H J = yn ® - ® V ° N ■ ( B1 ) 

Its dimension is obviously: 

JV JV 

<i j ]= e n( 2 ^ +i )= e n^®' (B2) 

T,3i = J i J2h = J i 

which can be interpreted as the dimension for the intertwiner space, but trading the group averaging integral over 
SU(2) for a localized evaluation on the identity g = I. The generating functional can be computed similarly. We 
thus define and compute the various generating functionals summing alternatively or simultaneously over the number 
of faces N and the total area J: 

dP[a,J] = Y, aN d°N[J}, (B3) 

JV 



N 



f N [t] = E* 2 '<M = ( E t 2j (2j + l) ] = [j^) , (B4) 



JeN yeN/2 

F°[a,t] = £ a»f N [t] = J2t 2J dP[a,J] = E^TtAw = ^.^L - ^ 

JVSN J N V > 

The lowest pole in t in terms of the variable u is easily identified, t c = 1 — and controls the asymptotics of d° [a, J}. 
Thus we get for large area J: 

S°[a,J] = log d° [a, J] ^ — 2 J log(l — ^fa) + . . . (B6) 

This reproduces the exact same leading order as in the intertwiner counting with the SU(2)-invariance. Thus this 
leading order, which gives the area-entropy law, is not related in any way to the requirement of SU(2)-gaugc invariance 
but simply to the straightforward sum over all possible spins on the N faces on the horizon. 
One can go further and provide exact formulas for the dimensions cfjy [J] : 

AM - ( 2N + S *\*J\- X>" ( 2N \)' ~ 1 ) ■ <■») 

From these expressions, it is straightforward to check that the log-correction to the entropy vanishes for large J. 
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2. Weakening the Closure Constraint: Quantum ^-Closure 

Let us now introduce the quantum model corresponding to a single closure constraint along the z axis. This amounts 
to drop the SU(2)-gauge invariance and replace it by the requirement of gauge invariance under U(l). This actually 
corresponds to the initial black hole entropy counting in Loop Quantum Gravity based on the identification of a U(l) 
Chern-Simons theory living of the black hole horizon Q . 

The number of states for fixed number of faces N and total area J is now derived through a group averaging over 
U(l): 

N 

4 [J] = dh IXj*W. (B8) 

where U(l) group elements are generated by the su(2) generator J z , more explicitly h = e*^ *. We further define the 
sum d}[a, J] = J2ngn oi N d^f[J] over the number of faces N . As before, we define the generating functionals: 

fj.it] = ^t w 4W = / * E t2j X j (h)\ , (B9) 

JGN M*)) \ 3 eN/2 J 

F 1 [a,t}=J2t 2J Ya N d 1 N [J}=J2a N fj f [t}=J2t 2J d 1 [a,J}. (BIO) 

,/ NGN NeN J 

Writing the integrals explicitly in terms of the class angle 4> of the group element, we get: 

,i m [*" d<f> ( v 2j S w(2j+l)A N [ 2 * d<i> 1 mi . , 

Mtl ~J ^ yj-' sin^ J ~ J 27(^-2icos0+ir' (BU) 

F*[a,t]= r~<f - 2tC °^\ l) . (B12) 
1 J J 2tt (t 2 -2tcos(j)+l) -a y ' 

The main difference with the SU(2)-invariant calculations (|39f and (|4Tj) is the measure term, where the sin 2 <j> factor 
now drops out. This integral can be computed explicitly from a, t ~ 0: 

F 1 [a, i] = 1 + " = 1 + a . (B13) 

v/(a-(l + i) 2 )(a-(l-t) 2 ) v/i 4 - 2(1 + a)i 2 + (1 - a) 2 

As before the lowest pole in i at fixed a €]0, 1[ is easily identified at t c = 1 — ^/a, which leads to the same asymptotics 
at leading order for the entropy: 

S l \a] = \ogd\f\a, J] = log ^o^dj^J] ~ -2 Jlog(l - y/a) + . . . (B14) 

AT J>>1 

One can go further and (easily) identify a first order differential equation satisfied by the generating functional F 1 [a, t]. 
Indeed, calling T = t 2 , we have: 

(T 2 - 2(1 + a)T + (1 - a) 2 ) St^ 1 + (T — (1 + a))(F — 1) = . (B15) 

This naturally induces a recursion relation on the coefficient of the series expansion F 1 = j T J d 1 [a, J] : 

VJ > 2, (1 - a) 2 (J + l)d x [a, J+ 1] - (1 + a)(2J+ l)^ 1 ^, J] + jd^a, J- 1] = 0. (B16) 

Inserting the ansatz d}\a, J] ~ A J /J a for large area J in this recursion relation, we get equations for the parameters 
A and a at leading and next-to-leading order in J: 

(1 - a) 2 A - 2(1 + a) + — = 0, (1 - ct)(1 - a) 2 yl - (1 + a) + a\ = . 



22 



These equations are very simple to solve and give: 



A=— J_, a=\, (B17) 
(1 - y/ay 2 



which gives the following asymptotics for the entropy: 

S 1 [a] = log 4 [a, J] j ~ i -2 J log(l - y/a) - 1 log J, (B18) 
with the expected factor — A for the log-correction. 

Appendix C: Bianchi's Semi-Quantized Model or Bohr-Sommerfeld Approximation 

As we see in Sec. Ill, the analysis for the classical polymer model for the density of states pn [A] corresponds to 
the first term in the expansion of the exact quantum result dim j [TV] for large area J at fixed number of faces N. An 
intermediate model was proposed in [l5j where the area of faces are quantized while the direction of their normal 
vector remains continuous and arbitrary on the 2-sphere. This half-quantized model, where the normal vectors do not 
become true quantum vectors but nevertheless have a discrete norm, can be seen as a Bohr-Sommerfeld approximation 
of the fully quantized model. 

This model is not particularly relevant to the present discussions, but we include its definition and a quick analysis 
for the sake of completeness. Hence we assume that areas of the individual faces are quantized Vj, = ami with 
m,i G N and a fixed area unit a (with dimensions restored) being of the order of the Planck area. But the directions 
of the individual vectors remain unrestricted. The total area is discrete in a-unit, now 2A = 2aJ. Then in this 
"Bohr-Sommerfeld" approximation the density of states is given as: 



n BS 

Pn 



m=4*E /(n^ 2 ^)^ 5 (E^- 2J ) 5(3) (E^)' ( C1 ) 



{mi} \i=l / \ k / \ k 



where a S- function on integers is the Kronecker 5. The the normalization is again chosen to enforce p BS [J\ = 1- 
Proceeding analogously to the derivation of Sec. Ill B I we write 

p B N s [J] =87rp ^-J £ J \J[m ia 2 e- am ^d 2 v}j e --£ fe ™*V a ^™^- 2J ). (C2) 



This leads to: 



{nii} \i=l 



P B n S [J] =**r£f ^K{a,ufe-^ , (C3) 



where the new kernel K(q, u) results from the integration over the normal vectors: 



oo „ 72tV • 2 

K(q,u) = Y m fl 2 e- f -e"" ra / UL e M ^ = h 2H^ a ( C 4) 

^— ' / 47r au 2(cosh(a(6 — in) ) — cos au) 

We haven't found a systematic way to perform the exact integration over the Lagrange multipliers q and u for all N 
and J so we are unable to provide an explicit closed analytic formula for the density of states in this model. We can 
nevertheless investigate the leading order effect of the area quantization in the N <C J — > oo regime where we keep N 
fix while considering large total areas J. 

Taking the TV-power of sm.au/au ensures that the dominant contribution to p^ s [Af] comes from ayfNu < 1, 
allowing to disregard the poles at u = q + 2Tm/a + ie for n > 1. In this regime we approximate the new kernel as: 

K(q, u) = 2 1 2 -$: + 0(a 4 ) « I(q, u) - a 2 /12. (C5) 
Then the BS density of state with the first non-trivial correction in the area unit a is given by: 

r?S\ n « 8^ ^ *L f l^ e ~2ia q j ( 1 N 



p N [J]~8nl —J—e ^ [ {u 2_ iq + u) 2 ) N i2(ti2_(g + ie )a)JV-i )■ 
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Performing the integration (and taking the limit ir/a — > oo) we find 

/ m3\ 

P % S [A = 2aJ] « p N [A = 2aJ] (l -—)+.... (C7) 

This expression obviously breaks down when N ~ J 2 / 3 , but we insist that this approximation holes for fixed N while 
J is sent to infinity. So this quantized area model obviously fits with the classical model at leading order, but its first 
correction in a (beyond the leading order in N <C J) unfortunately does not fit with the exact quantum expression 
(Sec. Ill), and thus this BS model will not be particularly relevant to our present work. 
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